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We study dynamical chiral symmetry breaking (DCSB) in an effective QED3 theory of d-wave 
high temperature cuprate superconductors under a uniform magnetic field. At zero temperature, 



'^ , the external magnetic field induces a mixed state by generating vortices in the condensate of 

O I charged holons. The growing magnetic field suppresses the superfluid density and thus reduces 



the gauge field mass which is opened via the Anderson-Higgs mechanism. By numerically solving 



^ . the Dyson-Schwinger gap equation, we show that the massless fermions acquires a dynamical gap 

m 

CN \ through DCSB mechanism when the magnetic field strength H is above a critical value He and the 

^^ ' fermion flavors A^ is below a critical value Nc- Further, it is found that both Nc and the dynamical 

m ■ 

^^ ■ fermion gap increase as the magnetic field H grows. It is expected that our result can be tested in 



phenomena in high temperature cuprate superconductors. 
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Analogous to Quantum Chromodynamics (QCD), Quantum Electrodynamics in (2+1) 
spacetime dimensions (QED3) has some interesting features, such as dynamical chiral sym- 
metry breaking (DCSB) in the massless fermion limit and confinement [lH9|. Therefore 
QED3 can be regarded as a useful model by which one can gain important insights into the 
aspects of QCD. In addition, QED3 and its variants also have proved to be the effective field 
theories for some two-dimensional condensed matter systems, including high temperature 



19|-|27|. Because of these features, QED3 has been 



superconductors |10l4l8l| and graphene 
extensively studied in recent years. 

DCSB occurs when the massless fermions appearing in a Lagrangian that respects chiral 
symmetry acquire a nonzero mass through nonperturbative effects at low energy. In a 
four-fermion interaction model, Nambu and Jona-Lasinio [28] first transplanted the fermion 
pairing idea in the BCS theory of superconductivity into particle physics to generate a 
fermion mass without assuming a Yukawa type coupling between fermion and Higgs scalar 
field (using Yukawa type coupling). In Refs. |3|, Ij] Appelquist et al. first studied DCSB in 
massless QED3 with N fermion fiavors at zero temperature by solving the Dyson- Schwinger 
equations (DSEs) for fermion self-energy in the lowest-order of 1/A^ expansion, and found 
that a fermion mass was dynamically generated when N is less than a critical number A^^ = 
32/7r^. Later, Nash further considered higher order corrections, and showed that the critical 
number of fermion fiavor still exists and Nc = ^ (5|. At finite temperature, Aitchison 
et al. used an approximate treatment of DSEs for the fermion self-energy in the 1/A^ 
expansion 29| and found that N^ is temperature-dependent and chiral symmetry is restored 
above a certain critical temperature. Further analogy with QCD at finite temperature 
and chemical potential has been explored 30l-l33J . The methods described there could be 
employed to extend the study of QED3 at nonzero temperature and chemical potential and 
explore whether and under which condition chiral symmetry restoration and deconfinement 
are coincident. Recently, the authors in Ref. 3^] studied the impact of a uniform magnetic 
field on explicit chiral symmetry breaking in QED3. 

Besides its interests in high energy physics, QED3 also turned out to be an effective 
low-energy theory of high temperature cuprate superconductors. In the past twenty years, 
it has been extensively applied to understand some unusual physics of high temperature 



superconductors [13|. In two early papers 



35, 



it was shown that the half-filled state on 



the Cu02 plane of undoped cuprates can be mapped to a (2+l)-dimensional continuum field 



theory composed of massless Dirac ferinions, called spinons, and U(l) gauge field, which is 
exactly QED3. In this effective theory, the physical fermion flavor is A^ = 2 < Nc-, so 
DCSB can take place at zero temperature. Ph ysic ally, DCSB corresponds to the formation 



of long-range antiferromagnetic (AFM) order [37|, |38|. At finite doping concentration, the 
simple QED3 of Dirac fermions is not sufficient to model the dynamics of doped holes. To 
overcome this difficulty, Lee and Wen developed an SU(2) gauge field theory to describe 
the underdoped cuprates [l3| . In the so-called staggered flux phase, two components of the 
SU(2) gauge field are gapped and thus can be neglected, leaving only one massless component 



of gauge field [13| . However, in underdoped cuprates, there appear additional scalar boson 
excitations, called holons, which describes the dynamics of doped holes. Therefore, the 
effective continuum theory of underdoped cuprates is an extended QED3 theory of massless 



spinons and scalar holons 37|, |38|]. The spinons and holons both interact with the gauge 

field strongly, but there is no direct coupling between them. This theory can exhibit a 

ot of interesting features, and has been investigated extensively in the past several years 
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Within this effective theory, superconductivity is achieved by the condensation of holons 
at low temperature. In the superconducting (SC) state, the gauge field acquires a finite 



mass gap through the famous Anderson- Higgs (AH) mechanism. As shown in Ref. [40||, this 
gauge field mass weakens the gauge interaction and thus rapidly suppresses DCSB. If we 
use DCSB to represent AFM order and use gauge field mass to represent SC order, then 
the competition between AFM order and SC order can be understood by the competition 
between DCSB and AH mechanism [40(]. 

Now suppose we introduce an external magnetic field into the SC state, then a mixed 
state, or called vortex state, emerges since high temperature superconductors are known to 
be extreme type-II superconductors. The most important influence of external magnetic 
field is to suppress the superfiuid density of the SC state. From a field-theoretic point 
of view, this implies the suppression of the gauge field mass, which is proportional to the 
superfiuid density. Since the gauge field mass can have important influence on DCSB, it is 
interesting to study the fate of DCSB when the strength of external magnetic field is varied. 
The main purpose of this paper is to study the impact of the external magnetic field on the 
critical fermion flavor A^^- 

In this paper, we first employ the DSEs method to construct the fermion gap equation 



at zero temperature after including the effect of an external uniform magnetic field. Then 
we solve the gap equation to study the influence of external magnetic field on DCSB in the 
mixed state of a d-wave cuprate superconductor. The main result of our paper is that both 
the critical fermion flavor Nc and the gap increase with growing magnetic field. 

It is well known that QED3 is super-renormalizable and has an intrinsic mass scale given 
by the coupling constant a = ^^. In Euclidean space, the Lagrangian of massless QED3 
with N fermion flavors reads 

£ = EV^.(^ + ie ^)^^ + -F% + -id, a,)\ (1) 

where the 4x1 spinor ipi represents the fermion field, i = 1,- ■ ■ ,N are the flavor indices, 
and ^ is the gauge parameter. Here, the massless Dirac fermions are not the electron-like 
excitations, but are actually the gapless spinous excited from the nodes of d-wave energy 
gap. These fermions have spin 1/2 but are electrically neutral. The gauge field a^ is also not 
the usual vector potential of electromagnetic field. Instead, it emerges as the consequence of 
strong correlation between electrons on the CUO2 planes of high-T^ superconductors. This 
effective field theory can be derived from the microscopic t-J model of high-Tc superconduc- 



tors with the help of slave-boson technique [13|. 

The DSEs for the fermion propagator G{p) at zero temperature reads 

^^^j^G{k)T,{p,k)D,M, (2) 

where Dpi, is the propagator of the a^ field and q = p — k. For simplicity, in this paper we 
adopt the rainbow approximation, i.e., we approximate the full vertex function ri,(p,k) with 
the bare vertex 7j, and we neglect the wave function renormalisation (this is a commonly 



emp 



19 



oyed approximation in the study of fermion propagator of QED3 at finite temperature 



41 



42|). Similar with the usual DSE study, in the following we will work in Landau 
gauge. 

At finite temperature T the DSEs for the fermion propagator is given by 

„2 00 „ rtr, 

G~\p) = Goip) - - E /T^^7'^G(fco,i^,/3)A^.(go,Q,/3)7^ (3) 

where 

Gip) = -■ ^ r^ (4) 

ipojo - vfPi'Ji - f AP272 -m{p^) 



is the ferinion propagator at finite temperature and m{p) represents the Dirac fermion 
mass. Here vp is the Fermi velocity, "^a is the gap velocity. Since the velocity anisotropy is 
unimportant to the critical behavior [l8|, we can set vp = V/\ = 1, and 



P= {Po,p),P = \p\,Po = (2m + 1 



TT 

k = {ko,k),K=\k\,ko = {2n + l)^, 

TT 



Q = (Qo, q), Q = |q| = |p - k|, go = 2(m - n)-. 
Under the above approximation the inverse of the fermion propagator can be written as 



(5) 



G (p) = i'y ■ p + m{p ). 



(6) 



Taking the trace of both sides of Eq. (3) yields a closed integral equation for m(p) 29|, |4l|, |42| : 

°° " d^k . , m(A;2) 



a 



H n=—cx 



(2vr) 



rA(go,Q,/3) 



^2 _j_ ^(J^2^ 



with 



where 



and 



A(go,(5,/3) = -Tr[7^A^^(go,<5,/3)7^ 



^^,v{qa,Q,l^ 


.X ^/.Sf^i^S 


' g2 + no(g,/3) + m2 


ow./?)^^; 


'^ 161n2 , 7r/3 ," 
TT 16 In 2 



(7) 



(8) 



(9) 



(10) 



When the photon acquires a mass ma via the Anderson-Higgs mechanism, the massless 
photon propagator is replaced by the massive photon propagator (9). So we have 

°° " cPk m{k^) 



m{p 



a 



E 



W/3 ^^^J {2ny [K' + m\k^) + wlW + MP. Q) + m'a\ ' 
The summation over infinite Matsubara frequencies can be done analytically: 

1 



(11) 



1 



= V 

2^K^ + m?{k^)nt^oo 

— . tanh 

2Ji\:2 + m2(A;2) 



1 



1 



iwn + JK^ + m^{k^) iwn- JK"^ + ni?{k 



-^K^ + m?{k^) 



(12) 



Finally we obtain the temperature- and momentum- dependent gap function 



"*P ' UNir'' J ^A'2 + ,„2(p) Q2 + n„(/!, Q)+my ^ ^> 

Here we note that in the limit of zero temperature /3 — )■ oo, we have 

no(g,/3)^^, (14) 



and the gap equation further simplifies into 



-^''^^^I''''jwBmQ'TWT<- '^«' 



In order to see how the gauge field acquire a mass m^, we now introduce the additional 
interaction term between gauge field and scalar boson field 

U = r (5o + mo)0 - ^r (5 + ^a)20 + ^02 + ^^4^ (ig) 

which is so-called Abelian Higgs model. The scalar field (p represents the bosonic holons, 
which has spin-0 and carry charge — e [l3|. This Lagrangian describes the motion of the 
charge degrees of freedom of electrons on the Cu02 planes of underdoped cuprate supercon- 
ductors [IS*]. Unfortunately, there is currently not a clear understanding of the dynamics of 



holons . In previous treatments of this model, the holon sector was simply neglected [37 - 
in the analysis of the interaction between Dirac fermions and gauge field. In this paper, we 
believe that the most prominent effect of this additional interaction is to generated a finite 
mass for the gauge field in the SC state. When m > 0, the system stays in the normal state 
and the vacuum expectation value of boson field (0) = 0, so the Lagrangian respects the 
local gauge symmetry. When m < 0, the system enters the SC state and the boson field 
develops a finite expectation value (0) ^ 0, then the local gauge symmetry is spontaneously 
broken. According to the Anderson-Higgs mechanism, the gauge field acquires a finite mass 
rria after absorbing the massless Goldstone boson generated due to spontaneous breaking of 
continuous gauge symmetry. The finite gauge field mass is able to characterize the achieve- 
ment of superconductivity. Actually, this mass is proportional to the superfluid density of 
the SC state, nia oc {(p). 

The external magnetic field turns the uniform SC state to a mixed state, where the 
superfluid density and therefore the gauge field mass decreases with growing field strength 



H . If we can get a relationship between gauge field mass iria and field strength H, we will 
be able to study the impact of magnetic field on DCSB. Since the holon djTiamics is poorly 
understood, we will use phenomenological arguments. 

As mentioned above, DCSB can describe the AFM order and the gauge boson mass 



can describe the SC order 
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40|. Based on these correspondences, we are able to give a 
qualitatively physical description of the competition between AFM order and SC order in 
the mixed state of high-Tc superconductors. In the absence of external magnetic field, the 
system stays in a uniform SC state at zero temperature, with a superfiuid density A^. As an 
external magnetic field is introduced, the phase coherence among holons is disputed and there 
appear single vortices inside which superconductivity is destroyed. In this mixed state, the 
superfiuid density becomes inhomogeneous and takes different value at different positions. 
However, after averaging over the vortices, it is still possible to obtain a homogeneous 
superfiuid density As{H), which is a function of magnetic field. On physical grounds, it is 
easy to suppose that As{H) is a decreasing function of field strength H. The gauge boson 
mass in the mixed state is proportional to the superfiuid density, so rria = aAs{H), where 
a is a dimensionless constant which depends on the doping concentration. Now we have 
a clear picture: as magnetic field H grows, the superfiuid density As{H) decreases, and 
hence the gauge boson mass rria decreases, which will raise the critical fermion flavor Nc 
and may finally induce DCSB in the mixed state. Once a finite fermion mass is dynamically 
generated, its magnitude m(p^) will be an increasing function of magnetic field H. In the 
following we will verify this qualitative argument by numerically solving the gap equations 
in the presence of an external magnetic field. 

To study the gap equation, the superfiuid density AJH) can be obtained by averaging 



over the vortices. In the study of mixed state, Volovik 43| at first proposed a semiclassical 



approach and showed that the density of states goes as yH at low temperatures, which has 



been observed by experiments 4^. A more complicated computation of superfiuid density 



As (.ff) within the semi-classical approximation has already been performed by Sharapov et 



al. |45|. The superfiuid stiffness is given by 

A'^{T,H) = T'^~Aii{T,H), (17) 

where Tij is the diamagnetic tensor and A„ represents the normal fluid density divided 
by the carrier mass. After taking some approximations and performing some complicated 



calculations, Sharapov et al. obtained 45 1 



^"'^'^^Tra:- 



When the superconducting gap Ad^^{H) <^ vH, we can neglect the A^^^ and A„(if) still 
has the behavior ~ ^/H. The superfiuid density then becomes 

As{H) = r- ^^Eh. (19) 

TCVa 

Here, r is the superfiuid density at zero temperature in the absence of magnetic field and 
the ratio between vp and f a appears as a coefiicient. When vp and f a multiply particle 
momenta, the anisotropy is irrelevant, but the anisotropy is very important in this expres- 



sion. In the following discussions, as in Ref. 45|, we simply set the ratio vp/v/^ = 30, and 



at zero temperature the superfiuid density is taken to be r = 1500K and the energy scale 



Eh[K] = 38KT~^/^wiJ[T]. Since the gauge boson m^ is proportional to superfiuid density, 
we use the ansatz that nia = aAs{H), where parameter a is a free variable. 

After obtaining the effective gauge boson propagator (15), we can now obtain the fermion 
gap equation in the presence of an external magnetic field 



The iteration algorithm can be employed to numerically solve this equation. In our numerical 
computations, we set a = 1 so that each quantity with mass dimension is re-scaled and 
becomes dimensionless. In the real numerical calculations, we adopt an ultraviolet cutoff 
A = 10^, which is large enough to ensure that the calculated results are stable with respect to 
A. We consider a = 0.07,0.1,0.15. For each parameter a, the relation between the critical 
fermion flavor Nc and the magnetic fleld H is shown in Fig. 1. The critical number of 
fermion flavor Nc separates the chiral symmetric phase (m = for A^ > Nc) from the DCSB 
phase {m > for A^ < A^^)- From Fig.l it can be seen that for flxed value of the parameter 
a, Nc increases as the magnetic fleld H grows. For a flxed magnetic fleld H, the smaller is 
the parameter a, the easier is it to generate an energy gap by the mechanism of DCSB. For 
the physical flavor N = 2, when the magnetic fleld H is large enough, a flnite mass gap is 
generated in the mixed state. It is obvious that the dynamical fermion mass will affect all 
observable physical quantities. It is expected that our result can be tested in phenomena in 
high temperature cuprate superconductors. 
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We now take the physical fermion flavor N = 2, and numerically solve the gap equation by 
iteration method for the specific case a = 0.07. The momentum dependence of the fermion 
gap is plotted in Fig. 2. From Fig. 2, it can be seen that the dynamically generated mass 
m(p^) is almost constant for small momenta, but decreases down to nearly zero at large 
momenta for a fixed magnetic field. This is easy to understand because the QED3 theory is 
asymptotically free. Another important feature is that fermion gap m{H) is an increasing 
function of the external magnetic field H. Apparently, the DCSB that is suppressed by finite 
gauge boson mass is restored when the magnetic field H becomes sufficiently large. When 

p2 

the fermion gap becomes larger with magnetic field, the thermal conductivity k oc ^a(-m i p^ ; 
where Fq is the disorder scattering rate, becomes smaller. There will be a thermal metal-to- 
insulator transition with growing magnetic field in the limit of m(H) ':$> Fq. These results 
qualitatively agree with the conclusion given by Jiang et al. in 46|, where the authors 
introduced a two-dimensional Coulomb interaction between Dirac fermions to generate an 
excitonic gap. 
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FIG. 1. The dependence of the critical number of fermion flavors N^ on the magnetic field for 
several values of a 



We note that QED3 in a uniform magnetic field was studied previously in several papers 



47H50|. A magnetic catalyst phenomenon was found there, which states that DCSB can 
be significantly catalyzed by an external magnetic field. In the present paper, this effect is 
not considered. This is mainly for two reasons. First, the gauge field in Eqs. (1) and (16) 
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FIG. 2. Scaled mass gap versus scaled monientum at various magnetic field for N = 2 and a = 0.07 



denoted by a^, is not the electromagnetic gauge field, but some new internal gauge field. In 
this sense the fermion does carry a gauge charge which is not the conventional electronic 
charge. Moreover, unlike the ordinary electrons that carry both spin and electric charge, 
the Dirac fermions appearing in our model are neutral in electric charge, so they are not 
directly coupled to the external magnetic field via the minimal coupling principle. Second, a 
crucial assumption of the magnetic catalyst phenomenon is that the external magnetic field 
induces discrete Landau levels of the Dirac fermions 



47 



49 



50|. However, it was found that 

51. 



there are no Landau levels in the mixed state of high temperature superconductor 

To summarize, there is a long history to solve the Dyson-Schwinger gap equation in 
QED3. It is interesting to apply it to some two-dimensional condensed matter systems. In 
this paper we apply the analysis to chiral symmetry breaking in high temperature cuprate 
superconductors under a uniform magnetic field. It is shown that chiral symmetry breaking 
takes place when the magnetic field is above a critical value He and the fermion fiavors 
is below a critical value Nc. Though the well-known magnetic catalysis does not occur in 
this model, the magnetic field leads to chiral symmetry breaking indirectly by suppressing 
the gauge boson mass. It is expected that our result can be tested in phenomena in high 
temperature cuprate superconductors. 
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